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Abstract 

We study nuclear structure effects on the transparency in high transverse 
momentum (p, 2p) and (e, e'p) reactions. We show that in the DWIA-eikonal 
approximation, even when correlations are included, one can get a factorized ex- 
pression for the transparency. This depends only on the average nucleon density 
p(r) and a correlation function. We develop a technique to include correlations 
in a Monte-Carlo Glauber type calculation. We compare calculations of T using 
the eikonal formalism and a continuous density, with a Monte Carlo method 
based on discrete nucleons. 

1 Introduction 

High p t (p, 2p) interactions in nuclei have gained renewed interest following the sug- 
gestion [1]] that fluctuations in nucleon size may show observable effects due to color 
screening. Of particular interest is the nuclear transparency, T, which is denned as the 
probability that a proton traversing a nucleus makes one and only one collision with a 
nucleon and that the emerging nucleons make no further interactions. It can apply to 
any interaction, but in this paper we consider the case where the basic interaction is ei- 
ther a proton-proton or an electron-proton interaction. The transparency is extracted 
from the experimentally obtained quasielastic (p, 2p) events using the factorization 
assumption |2|, |3|: 

d r -» d 

—<j p - A /Z = T J S(k,e,a tot )—a p ^ p (s,t)d 3 kde (1) 

This is the expression found in Eq. (2) of Ref. ||. In this factorized form o" p _ p is the 
(p, 2p) cross section evaluated at the appropriate Mandelstam variables, s and t. The 
variable s(k, e) is the on-shell center-off-mass energy which is a function of the missing 
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energy, e, and missing momentum, -k, measured from the kinematic reconstruction of 
the (p, 2p) events, and t is the measured 4-momentum transfer. Finally S(k, e, a to t) is a 
spectral function, which depends on the proton-nucleon interaction cross section, a tot , 
and on the spatial distribution of nucleons in the nucleus. It is normalized to unity. 

Using Eq. ([]]) T can be extracted from the experimental data by weighting each 
observed event with the reciprocal of ^cr p _ p (s, t), to remove the dependence of T on k 
and e. For an experiment designed to accept events of any momentum k and missing 
energy e there is no need to know S(k, e, <J to t), since / S(k, e, a to t)d 3 k de = 1. In practice, 
at incident momenta of 6-10 GeV which are much larger than the average k, corrections 
for apparatus acceptance as a function of k are not difficult to obtain, especially since 
5* falls rapidly with k. 

If there were no initial state or final state interactions of the protons traversing the 
nucleus, T would be unity. Thus T measures these initial and final state interactions 
which are sensitive to the total proton-nucleon cross section, a to t, which in turn will 
depend on the magnitude of the color screening. One of the important ideas in color 
transparency is that T will increase with the momentum transfer, t, observed in the 
(p, 2p) reaction. Thus we examine how to determine T(t) in a way that will depend least 
on theoretical calculations of its magnitude. This is also of value since it bypasses the 
strong sensitivity of T to the nucleon density near the nuclear edge. (In general only a 
peripheral reaction will have a single scattering, so that edge effects are enhanced.) We 
examine the special case of the t dependence of the nuclear cross section as measured 
in a (p, 2p) reaction where the exiting protons are observed at the angles 9\ and 82- 

In this paper we show that a factorized expression for T, Eq. (0), can be derived 
from the Distorted Wave Impulse Approximation. By using the closure approxima- 
tion, as employed by Lee and Miller |4[], the spectral function S becomes independent 
of e (we denote this simplified function by J in this work). Even with this additional 
approximation, the expression for the cross section is rather complicated, see Ref. 0. 
However, we show that the transparency depends only the nuclear proton radial dis- 
tribution p(r). The density p can be obtained from elastic e — A scattering without 
recourse to knowledge of the individual single particle wave functions. After using clo- 
sure the transparency T turns out to be related to the classical probability that there 
are no initial or final state scatters accompanying the (p, 2p) high p t interaction on a 
single nucleon. We shall see later how to take into account both the momentum and 
missing energy taken off by the undetected nucleons measured in a fully reconstructed 
(p, 2p) experiment. 

We also show how the effect of nucleon correlations, which affect the calculation 
of T, can be included in a previously developed Monte-Carlo formalism for Glauber 
scattering, and compare various model assumptions with the data. 

We also discuss the transparency in the (e, e'p) reaction, for which both theoretical 
H and experimental studies |]] have been performed. The formalism developed by 
Benhar et al. closely parallels the DWIA eikonal formalism used by Lee and Miller. 
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2 Transparency and the Distorted Wave Impulse 
Approximation 



Let us first review the expected cross section on the basis of the plane wave impulse 
approximation (PWIA). This cross section can be obtained by using the closure approx- 
imation to sum over final states and using on-shell values of the Mandelstam variables, 
s and t, for the quasi-free interaction 0. The result of this quantum-mechanical ap- 
proximation is to obtain a simple limit: The cross section, for an incident proton of 
energy, £o, on a nucleon with wave function ip becomes the probability, P{k), of find- 
ing that nucleon in the nucleus with momentum k, times the (p, 2p) cross section at 
the on-shell values of s(Eq, k, e) and the momentum transfer, t. While t is directly 
obtainable from the measured momenta, k and e refer to the momentum and binding 
energies in the initial state which are not directly measurable. 

For interaction with a single nucleon of independent particle wave function t/>, the 
result is well known ||: 

^p±/Z = [p(k) ^-* [s } k > € > t)] <Pk. (2) 

(JjL «/ (Jju 

In particular P(k) = \F(k)\ 2 where F(k) is the Fourier transform of the nucleon 
wave function, i[)(r), 

F(k) = J e* r >(f)d 3 r. (3) 

It is useful to recall that, in Eq. (|3[), the wave function, ip(r), must be normalized 
to unity in order that / P(k)d 3 k give the proper total probability, i.e., J P(k)d 3 k = 1. 

We next examine the cross section in the distorted wave impulse approximation 
(DWIA). In this approximation the plane waves for the incoming (i = 0) and outgoing 
(i = 1, 2) particles are replaced by distorted waves of the form e 1 ^'^ Di(r, p(r), a to t,Pi), 
where Di is a factor describing the distortion of the plane wave. It depends on the 
total p-nucleon cross section a t ot, the nuclear density p and the momentum of the ith 
proton pi. 

For the case of a single filled orbital we find that F(k) in Eq. @ is replaced by 
F'(k), 

F'(k) = J e^DoD^i/jirj^r. (4) 

Note that since in general the D's are not rotationally invariant the momentum dis- 
tribution no longer depends only on the magnitude of k but depends on both the 
transverse momentum, k t , and longitudinal momentum, k z . 

In order for P'{k) = \F'(k)\ 2 to be a proper effective momentum distribution which 
includes the initial and final state interactions, F'(k) must be normalized so that 
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Figure 1: Geometry of the scattering process. The vectors b{ are the impact param- 
eters for each of the three paths (0 labels the incoming path, 1 and 2 label the two 
outgoing paths) for a collision taking place at a nucleon at r. z , and equivalently z\ 
and z 2 , are the z coordinates of r in the b, z coordinate systems. Q\ and 9 2 are the 
angles of the outgoing protons relative to the incoming proton direction. 



/ P'(k)d 3 k is unity. The normalization factor needed in Eq. (f|) is defined as VT, 
anticipating the result displayed in Eq. fl7|), and is given by: 



D D 1 D 2 ^(r)*D D 1 D 2 i;(r)d 3 r, 



so that 



P'(k,a) 

We then find, in analogy with (0 



J e i ^D D 1 D 2 ip(f)d 3 7 



T 



(5) 



(6) 



da, 



p-A 



dt 



/Z = T J P'(k) 



da, 



p_p L s(fc, e, t)} 3 ^ 



dt 



(7) 



The p-p cross section can be obtained from p-p data where s is determined from 
k and e, which in their turn are obtained by kinematically reconstructing each event 
under the assumption that all the accepted events are pure (p, 2p) events with no initial 
and final state rescatterings. 

In order to obtain the DWIA expression for T we follow Lee and Miller and use the 
eikonal form of D to evaluate T. With the path lengths, z, defined in Fig. |T] (z along 
the direction of the incoming particle, z\ and z 2 along the trajectories of the outgoing 
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particle) we have: 



Di 
D 2 



exp 
exp 
exp 



—p(z,b )dz 

-oo Z 



^p(z' 'M)dz' 



&tot 



I 22 2 

Inserting these expressions in Eq. (ffi) we obtain: 



p(z", b 2 )dz h 



(8) 



Jl x a totP (r)dz 0e - £» <x tot p(r)d 2le - /« = j p^p^^ ( 9 ) 



The functions P,, that are the squares of the functions -Dj in (^), can be interpreted as 
the probability that there are no interactions in the leg i of the reaction. 

In general we have several (partially occupied) orbits in a shell-model nucleus. The 
result found here remains true, however. Let us consider the example of two completely 
filled levels, s and p, with equal occupation numbers. We define the "normalized to 
unity" effective momentum distribution for state i (i = s,p) by 



fi(k) 



in. 



with the normalization factor determined from: 



Ti = [(D D 1 D 2 ip i (?))*[D D 1 D 2 ,l; i (f?))d 3 



We then have the useful form: 

a p - A /Z = (T s + T p ) J 



fsT s ~\~ fpTp 



T s + T v 



x 0p-p(s, t)d 3 k. 



(10) 



(12) 



Note that T s + T v is identical with T since 



T s + T p = J (D DiD 2 ) 2 [ip s (ryil; s {r)+'4)p(ryip p (r))d 3 r = J (D D 1 D 2 ) 2 p(r)d 3 r. (13) 

This result can easily be extended to more independent particle states and is the 
basis for our result that only p(r) is needed to calculate T. Thus Eq. ( |T3"D is identical 
with Eq. (|). 

Note also the useful relation for the momentum distribution, resulting from our 
normalization of /: 



J(k,a)d 3 k 



fsT s + fpTp 
T s + T p 



d 3 k = 1. 



(14) 
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The final form can thus be written: 

j° v -aIZ = T J J(k, cr)j t a p . p (s, t)d 3 k. (15) 

Thus, as in the PWIA, the DWIA also leads to a cross section that can be factored into 
an effective normalized momentum distribution, J(k), and a k independent quantity, 
T, called the transparency. This result relies on the use of a closure approximation, so 
that there are no interference terms in Eq. (0). Let us repeat this equation (|9]) once 
more, since it is one of the important results of this paper: T — J PoP\P2p{r)d 3 r. Thus 
we recognize that the normalization factor in Eq. (^) is in fact just the transparency 
in Eq. (||)! This equation is exactly the expression for T used by Farrar et al. in 
their calculations. 

Eq. (0) can be understood as the probability that there is no nuclear interaction of 
the incoming or outgoing protons along the paths Zq , z\ and z 2 determined by r, 9i , and 
6 2 . However it does not include "nuclear correlation" effects that modify p(r) in Eq. (0) 
in the neighborhood of the struck nucleon. These corrections were calculated in Ref. || 
and are discussed and included by Lee and Miller M. We have also calculated this 
probability using a Monte Carlo method which automatically || excludes the nucleon 
participating in the high p t (p, 2p) interaction from the absorptive path. 

Equations (§) and fllS]) allow T to be determined directly from the proton charge 
density, p(r), which can be obtained from elastic e- nucleus scattering rather than 
requiring a detailed knowledge of the ground-state wave functions for all the nucleons 
of a complicated nucleus, which is the procedure followed in Ref. Q. 

The correlations of the struck nucleon with its surroundings make it less plausible 
to find another particle very close to this nucleon. This can be most simply be taken 
into account by replacing the one-body density p in all three equations (0) by the 
probability to find a particle at position along one of the three legs, if there is one 
at the beginning (at position r), 

P2(r, r zh )/p(r) = p(r zb )C(f, r zb ). (16) 

Here we have defined a correlation function C. This entity is usually approximated by 
the result of a nuclear matter calculation. As is argued in the appendix, however, this 
correlation function should also take into account the (A — 1)/A corrections arising 
from the fact that one should not include the struck nucleon among the absorbing 
material. Thus, for no correlations, we find C = (A — 1)/A, which for light nuclei gives 
a considerable correction to the transparency. 



3 The "Glauber" Calculation of T 

The Glauber model has proved to be quite successful in understanding p-A and heavy 
ion interactions at high energies (> 5 GeV). The basic quantities such a calculation 
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Figure 2: The geometrical cylinder where no particles are allowed for a transparent 
reaction. For the calculation of a n we draw a cylinder through the whole nucleus 
and count the number of particles inside. Here we have illustrated one of the legs in 
the (p, 2p) scattering, centered on the particle where the hard collision is made. We 
construct a similar cylinder for each of the three legs. 



are the so-called Glauber coefficients, a n . The coefficient a n gives the probability that 
a proton makes n collisions in traversing a nucleus, while remaining on a straight path 
(this assumption is similar to the eikonal approximation). One approximate way for 
obtaining the Glauber coefficients is based on a Monte-Carlo generations of particles 
in a nucleus. One first generates the positions of nucleons in a nucleus randomly, 
according to the nuclear density p(r). One then follows the path of the incoming 
proton through the nucleus, and counts the number of collisions. The collisions are 
treated geometrically, and are assumed to take place whenever the impact parameter 
for the incoming proton and any target nucleon is less than b a obtained from the proton- 
nucleon cross section, with <j tot = nb^. The number of collisions is thus equal to the 
number of centers of particles in a cylinder of radius b a around the path of the projectile, 
see Fig. Here we have taken into account the fact that at very high energies the path 
of the incoming particle through the nucleus can be approximated by a straight line. 
Such a probabilistic calculation of the a n is known to be equivalent to the quantum 
mechanical Glauber model in the case of factorized ground state densities |{J and has 
been used over the last decade to determine the a n which are used in models that study 
multiplicities [|H], |TT|, E t distributions [|1^, |TB|, leading hadron distributions [jl4[ and 



many other measured quantities in p-A and heavy ion collisions at high energies. 

For our transparency calculations we need determine the probability that if the 
projectile nucleon makes a large-angle (p, 2p) scatter, no other soft collisions take place. 
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Figure 3: The transparency calculated as a function of A for the Glauber Monte Carlo 
model. The three curves represent three different laboratory angles for the pp scatter. 

We sum over all particles in the nucleus, and assume that the hard-collision takes place 
at one particle at a time. We then look at the incoming path, and also follow the paths 
of the two final state protons, and look whether a geometric scattering takes place 
along each path. We allow the nucleons to scatter into non-zero angles as determined 
by the (p, 2p) kinematics. The calculation is similar but not identical with the eikonal- 
DWIA calculation discussed in the previous section. The expression for T resembles 
Eq. (^) with the exponential factor replaced by the probability that there are no 
other interactions except the hard pp interaction along the appropriate path. 

In figure |3| we show the dependence of the transparency on the laboratory angle of 
the pp scatter. The kinematics are usually chosen so that the cm scattering angle is 
90°, we see that the transparency increases trivially as the energy increases. We have 
used Utot = 37 mb, which we consider to be a reasonable value for the high-energy pp 
cross-section in nuclei. This value will be used in all calculations in this paper. 

One basic difference between this Glauber calculation and the DWIA calculation 
is that the Monte Carlo treats the nucleus as composed of confined nucleons of finite 
volume while the DWIA treats a continuous matter distribution. We believe that at 
10 GeV, where the wavelength of the projectile is smaller than the nucleon diameter, 
one cannot neglect the "lumpiness" of the nuclear medium. In contrast to the eikonal 
calculation discussed in the previous section, in the Glauber calculation the target 
nucleon participating in the hard (p, 2p) scatter is by definition not counted in the 
absorptive path. This does not mean that nuclear correlations do not have an effect on 
the result of the Glauber calculation, but the effect can be different. Let us investigate 
the generation of correlated nuclei using our Glauber-Monte Carlo algorithm. 
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3.1 Monte Carlo generation of particles 



As in the Eikonal-DWIA calculations we wish to use a nuclear matter calculation as a 
guide to construct a correlated sampling of a finite nucleus. There are several reasons 
why we cannot use a correlation function directly as input to the calculation. The most 
important reason is the fluid nature of the nuclear many-body system, which induces 
many-particle correlations. 



3.1.1 Nuclear matter 

Let us first look at the case of nuclear matter. If we consider particles with only 
two-body interaction, one expects that the probability of finding the particles can be 
written as a product of two-body correlation factors, 

p(n,...) = l[pl[g(\f l -f j \) (17) 

i i<j 

where g is a two-body correlation factor. This function should approach one for large 
separations. The point is that, for nuclear matter, the correlation length - the distance 
in which g(r) approaches one - is comparable to the internucleon spacing. This means 
that in general we will be able to find triples of particles (r*j, r*,-, r^) in the Monte-Carlo 
sample such that at least two of the three g's {gifi — fj), g(f\ — fk) and g(fj — ff.)) 
differ significantly from one. This leads to effective (induced) three-body correlations. 
If we now calculate the two-particle distribution, the probability to find one particle 
at fx and another at r 2 , 

P2(r 1 ,f 2 )=p 2 C(\f 1 -f 2 \), (18) 

we find that the function correlation function C is in general different from g. (Note 
that in the low-density (gas) approximation C = g.) This can be seen in Figs. |] and 
H] for two choices of the function g. We have used a Monte-Carlo algorithm where we 
generate N particles in a square box with periodic boundary conditions. 
In Fig. H we show the effect of a hard-core interaction, 



»m = : z ■ da) 




As one can see the resulting curves do not have very much in common with the results 
of more standard nuclear matter calculations (see Fig. |). The hard core reduces the 
available positions in configuration space so that for the large core radii each particle 
is surrounded by a number of particles that almost touch. This leads to the strong 
peak in the correlation function for r = R. In Fig. |5| we use a Gaussian g, 

g(r) = l-exp(-(r/R) 2 ). (20) 

This function results in a C that - for R = .92 fm - is very similar to the realis- 
tic correlation function used by Benhar et al. || in their calculation for the (e, e'p) 
transparency. This is what we use in our further calculations. 
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Figure 4: The effect of hard-core correlations in "nuclear matter", on the correlation 
function C, for various hard-core radii. The solid line is hard core radius. In sequence 
the other three lines have a hard-core diameter of .46, .92 and 1.38 fm. 



1.2 




r (fm) 

Figure 5: The effect of a Gaussian correlation in "nuclear matter", on the correlation 
function C. The dashed line is the input Gaussian, the solid line the correlation 
function used by Benhar et al. The black circles represent the outcome of a Monte- 
Carlo calculation. 
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3.1.2 Finite nuclei 



Whereas for infinite nuclear matter we can have an arbitrary, but translationally in- 
variant, correlation function, the correlation function for finite nuclei must satisfy a 
normalization condition. This follows from the equation 

r/ :; r 2 p 2 (n.^ ( A \),>{r 2 ). (21) 
which after the substitution 

p 2 (n, r 2 ) = p(ri)C(fi, f 2 p(f 2 ) (22) 

becomes 

J d 3 r 2 C(r u f 2 )p(f 2 )/(A - 1) = 1. (23) 

We obviously also require that the one-body distribution coming from the Monte-Carlo 
generation takes on the form p(r). 

We wish to develop an approach to finite nuclei that is close in spirit to the nuclear 
matter calculation discussed in the previous subsection. Therefor we use an approach 
based on two-body correlation factors. We use a iV-particle distribution p of the form 

N N 

PN {n, . . . , f N ) = ((n) HgQfi - r,-|). (24) 

i=l i<j 

M is a normalization factor, that could be absorbed in the effective single particle 
densities (. The function ( is chosen such that p{f) takes a given form (usually Wood- 
Saxon), whereas g is taken to be identical to the nuclear matter result. For a dilute 
gas one find that ( = p and g = C . Unfortunately a nucleus is more like a fluid than 
a gas, so that we have to solve for (. 

In this work we require that for g given by Eq. (^) we obtain a Wood-Saxon 
density, 

P( r ) oc — ^77^, (25) 

1 + exp([r - R \/cr) 

with a = .545 fm and Rq = 1.14v4 1//3 fm. It is rather complicated to solve the many- 
body problem for (. We found it convenient to use a parametrized form and look for 
a set of parameters that yields a result that closely resembles a Wood-Saxon. A form 
that seems to be sufficient for this goal is 

C(r) a (l + exp([r-^]/a'))(l + «r)' (26) 

The function £(r) is a parent distribution, from which we chose particles in the Monte- 
Carlo calculation. If we include the correlation factors g (by a rejection technique) this 
leads to a single particle distribution of the desired form p(r). 

To illustrate the accuracy obtained by using ( p6| ) to arrive at the Wood-Saxon 
distribution fl25), we show a few of the results of our calculations in Fig. || The 



relevant parameters for a few nuclei are given in Table [TJ 
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Figure 6: A comparison of the result of the density for the correlated Monte- Carlo 
result to a Wood-Saxon, for four different nuclei. For each pair of graphs the upper 
one shows the value of the single-particle distribution in a set of bins. The triangles 
are the Monte-Carlo results, and the drawn line is the exact result for a Wood-Saxon 
distribution. The lower graph shows the difference, A, between the two sets of results. 
For parameters see Table [1] 
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Table 1: The parameters of (, Eq. |26] for a few nuclei, as determined from a Monte- 
Carlo calculation. 

A R (fm) R' Q (fm) a' (fm) a (1/fm) 



7 


2.18 


2.18 


.545 


0.10 


12 


2.61 


2.55 


.545 


0.13 


16 


2.87 


2.75 


.545 


0.16 


40 


3.90 


3.70 


.545 


0.19 


100 


5.29 


4.90 


.540 


0.19 


208 


6.75 


6.17 


.530 


0.19 



4 Model Dependence of T 

To elucidate the sensitivity of T to the shape of the nuclear surface we show in Fig. |^ 
the dependence of N(b)/b, the number of collisions at impact parameter, b, versus b. 
Superimposed on the same figure is the plot for the number of pure (p, 2p) events, i.e., 
those without initial or final state scattering, T(b) (multiplied by 10). These curves 
illustrate how the tail of the N(b) distribution determines both the shape and position 
of T(b). Nost of the transparent events take place where the nuclear density is about 
1 /4 of the maximum density. 

In Fig. H we compare the results of various calculational schemes discussed in this 
paper. We have calculated the transparency in (p, 2p) reactions for the methods dis- 
cussed in the previous sections, all using a Wood-Saxon density of the form 

" W K l+expKr-flo)/,] <27) 

where Rq = 1.14A 1 / 3 and a = .545 fm. As stated before the total pp cross section was 
(in all cases) taken to be a rmta t = 37 mb. We have calculated the transparency for 
the DWIA without correlations, for the correlated DWIA-eikonal approximation, with 
the correlation function used by Benhar et al. - the one used by Lee and Miller gives 
almost identical results -, the Glauber Monte Carlo calculation without correlations 
and this last method with correlations included. The smallest transparency is found 
for the uncorrelated DWIA calculation, comparable to the work of Farrar et al. If we 
add correlations to the eikonal calculations, in the form of a nuclear matter correlation 
function, the transparency increases substantially. Similar results are found for the 
Glauber calculation, where the results are consistently higher than those from the 
eikonal calculation. The data of Carroll et al. |2[] are superimposed on the family of 
calculations in Fig |8]. 

To ease the comparison for large A we have also tabulated the transparencies in 
Table [2| From this table we see that the correlated eikonal is, even for high A, 30% 
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0.010 




b(fm) 

Figure 7: Impact parameter dependence of collisions in 64 Cu. Filled circles denote 
the number of all scatters vs. the impact parameter, b. Open circles give the number 
of transparent collisions multiplied by ten. Both figures use a Woods-Saxon density. 
The upper figure has very small diffuseness, o = .005 fm, the lower has a realistic 
diffuseness, a = .545 fm. 



Table 2: Transparencies a calculated for Fig. [8[ The column labeled DWIA gives 
the Farrar-like results, DWIA cor gives the correlated eikonal-DWIA results with a 
Benhar type correlation function. MC labels the uncorrelated Glauber Monte Carlo 
calculations, and MC cor gives the correlated result. 



A 


Tdwia 


^DWIA cor 


Tmc 


Tmc cor 


7 


0.342 


0.446 


0.419 


0.512 


12 


0.242 


0.327 


0.291 


0.365 


16 


0.197 


0.270 


0.233 


0.292 


40 


0.0965 


0.134 


0.108 


0.136 


56 


0.0732 


0.101 


0.0807 


0.100 


100 


0.0451 


0.0611 


0.0476 


0.0575 


208 


0.0243 


0.0319 


0.0254 


0.0298 
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Figure 8: The transparency for the (p, 2p) reaction calculated using various methods. 
The angles 9 have been chosen appropriate to a geometry corresponding to 90° scat- 
tering in the cm frame for Ei a b = 10 GeV. The 6, 10, and 12 Gev data are from 
Ref. §. 

larger than the uncorrelated eikonal calculation. The effect of correlations on the 
Glauber Monte-Carlo calculation depends somewhat stronger on A, and is 22% for 
A = 7 and 17% for A = 208. 

We have found an error in the calculation of Ref. [H. The effect of including 
correlations increases the transparency by 60% and not by 100% as previously reported. 
As can be seen from Table [| the Monte Carlo method gives somewhat larger values of 
T than the eikonal approximation for a similar correlation function. 

Since our work was completed an important calculation of T for the (e,e'p) inter- 
action has appeared in print ||. Benhar et al. have derived an expression for the 
transparency which carefully incorporates the best present knowledge of the nuclear 
correlations. The expression obtained in their work is very similar to the one derived 
in Sec. 2. We have performed similar calculations, corresponding to removing any of 
the particles in the nucleus in a given direction. For the eikonal calculation this leads 
to an expression 

T = J d 3 rPi(r)p(r), (28) 

and there is a similar simplification for the transparency in the Glauber scheme. In 
fig. [5] we compare the same calculational schemes as before. For reference we also show 
the the (p, 2p) results. The difference in curvature as a function of A can be explained 
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Figure 9: The transparency for the (e, e'p) reaction calculated using various methods. 
For ease of comparison we have also included the (p, 2p) results. 



as follows: imagine the nucleus as a large sphere. If a proton comes in parallel to the 
line connecting the poles, transparent events will only take place near the equator (on 
the surface). The fraction of transparent events is thus proportional to R/R 3 oc A~ 2 ^ 3 . 
For the (e, e'p) reaction where the momentum transfer is parallel to the same axis the 
whole southern hemisphere will contribute, leading to a transparency proportional to 
R 2 /R 3 <xA~^ 3 . 

It is in our geometrical approach very natural to exclude the struck particle from 
rescattering the projectile. In the original form of the eikonal model this is not done. 
One should replace the one-body density p in the exponential by the ratio of the two- 
body and one-body density functions. Lee and and Miller have modified their DWIA 
calculation to take into account this effect approximately. This is done by multiplying 
p(r) in Eq. @ by a nuclear matter correlation function C(z — Z{). The correlation 
factor reduces the path length in the vicinity of the point where the proton undergoing 
the (p, 2p) interaction is situated. They find that such inclusion increases the value of 
T by about a factor of 1.6 for Carbon. Thus their result supports the remark M that 
the Farrar et al. calculation neglects this important effect. While the increase in T 
due to the correlations is similar in both our calculations, the magnitudes may not be 
the same since the effective correlation functions differ in the DWIA and probabilistic 
approaches. Lee and Miller use a nuclear correlation function obtained from an earlier 
paper |T5| , which differs from the nuclear correlation employed by ref . p| . 

Actually our results labeled "Lee and Miller" differ from those in their paper (H) 
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Table 3: Calculated Values of the transparency T for 12 C in the (p, 2p) and (e, e'p) 
reactions. The kinematics used are appropriate to 90° scattering in the cm frame at 
a lab energy of 10 GeV. Lee and Miller denotes the correlation function employed in 
Ref. gj, Benhar that in Ref. §. 

Correlation 



Method 



p(r) 



T, 



T 

(p,2p) ± (e,e'p) 



Eikonal 


Woods-Saxon 




0.24 


0.60 


Eikonal 


Woods-Saxon 


Lee and Miller 


0.32 


0.67 


Eikonal 


Woods-Saxon 


Benhar et al. 


0.33 


0.68 


Glauber 


Woods-Saxon 




0.29 


0.63 


Glauber 


Woods-Saxon 


This work 


0.44 


0.72 


Eikonal 


Gaussian 




0.17 


0.51 


Eikonal 


Gaussian 


Lee and Miller 


0.24 


0.62 


Eikonal 


Gaussian 


Benhar et al. 


0.25 


0.63 


Glauber 


Gaussian 




0.27 


0.54 


Eikonal 


(e, e') 




0.18 


0.53 


Eikonal 


(e, e') 


Lee and Miller 


0.26 


0.62 


Eikonal 


(e, e>) 


Benhar et al. 


0.27 


0.63 


Glauber 


(e, e>) 




0.28 


0.55 



for another reason: The definitions of T are not identical. Lee and Miller define T as 
the ratio of the theoretical DWIA cross section to the theoretical PWIA cross section 
which is not the same as the expression used by the experimentalists (Eq. (|15|) ) to 
extract T from the data. As we can see from Eq. (0) and ( |T5D their ratio, Tim, is given 
by 



LM 



TJ J(k,a) d(7p ^ s ' t) d 3 k 



(29) 



Since J(k, a) is not equal to P(k), the Lee and Miller definition is not that used in the 
experimental analysis. 

The different predictions for Carbon are summarized in Table [| We have used a 
Wood-Saxon density as well as a Gaussian density, obtained from using a completely 
filled set of s and p orbitals with a length parameter 6 = 0, similar to the one employed 
by Lee and Miller. We have also used an experimental charge-density obtained from 
elastic electron scattering in these calculations |i~6|| . The results with this last density 



are close to the one obtained with the harmonic-oscillator. As can be seen for the (p, 2p) 
reaction the result of the correlated eikonal and the Monte-Carlo Glauber calculation do 
not differ very much. The effect of correlations on the Monte-Carlo calculations leads 
to a result 33% larger than the correlated eikonal calculation. This result seems to be 
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mostly due to finite A corrections, which are missing in the eikonal calculation, and are 
important for this light nucleus. It is important to note that the effect of correlations 
in the Monte-Carlo calculations is to significantly enhance the transparency. This 
shows that one should take into account two effects. First one should not count the 
participants of the hard pp scatter among the absorbing volume. Secondly, the effect 
of correlations on the local environment of the struck particle is ver important as well. 

All of these effects, the sensitivity of T to the shape of the tail of p(r), to the 
nuclear "granularity", or to the effect of "nuclear correlations" in the eikonal model 
serve to emphasize that the magnitude of T may be a poor quantity to use for color 
transparency studies. 



5 Extracting T from the Data 

We note from the defining equation, Eq. QT5|), that the cross section depends on three 
quantities, S(k), T, and a(s, t). If our goal is to search for the presence of color screen- 
ing, the t dependence of T, and not its magnitude, is the crucial quantity. Also, S(k) 
has intrinsic interest, but not necessarily to the extraction of T(t). 

The experimental process for extracting T from the data has been given in Ref. 
and also in the thesis of Guang Yin Fang [|Tj]] . The experimentalists take the factorized 
form (Eq. fll5D ) of the nuclear cross section as their basic assumption which we see is 
the same as the DWIA result. To extract T each observed event can be weighted with 
W s i = 4cr(0, 0)/4cr(A;, e) = ^a PjP / -^a(k, e) to correct for the cross section dependence 
on k and e, and is also divided by Z. (Ralston and Pire |TB|| have pointed out the need 



to consider variations of the (p, 2p) cross section from the simple s~ 10 behavior.) 

If there were no acceptance biases, so that one could use the fact that the integral 
over all k is unity, T would be obtained directly. 

Since, strictly speaking, T = T^i,^), and since events with different t might 
appear at different angles, one can consider adding to future analyses the theoretical 
weighting factor, = 1/T(9i,9 2 ). Otherwise a 9 dependence on t might appear as 
a spurious dependence of the transparency on t. The rise in T resulting in neglecting 
this effect can be read from Fig. ^ for the special case of an equal angle configuration. 

Although we have emphasized that the magnitude of T is model dependent, the 
dependence of T on the exiting angles would not be expected to be. (The weight, Wi t2 , 
is easily calculated from our available Monte Carlo code.) Thus we conclude that the 
data analysis can be refined by using this angular weighting factor for more precise 
determination of the average T in the experiment. 

It might seem on a casual analysis that one could use the transverse missing mo- 
mentum spectrum observed from the data to predict the longitudinal momentum dis- 
tribution, since the transverse momentum distribution is not sensitive to the strong s 
dependence of the cross section. However whenever initial and final state interactions 
depend on the projectile axis, the transverse and longitudinal momentum distributions 
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in J(k, a) are not identical. A nice internal check of the data would be to unweight 
the k z distribution with the cross section to see how much it differs from the trans- 
verse momentum distribution. This is important if the extraction of T depends on a 
knowledge of these distributions in order to correct for apparatus acceptance as well 
as momentum cuts that might be applied to the data. 

6 Discussion 

It is important to realize that the DWIA approximation results in a factorized cross 
section which involves an effective momentum distribution, J(k,a). The binding en- 
ergy e has disappeared from this expression because of the sum over all final states in 
the closure approximation. If the experiment were to detect all (p, 2p) events indepen- 
dent of k and e, this would not matter because that sum is unity. But S(k, e) is needed 
to carry out the Monte Carlo calculations used to find the true apparatus acceptance. 
This is just another way of saying that each event % must be weighted by a Wi(k,e). 
Since protons of the same k can be ejected from states with different binding energies 
(and there is really a spectrum of excitation energies for each initial state), it is clear 
that the S(k, a) appearing in the factorized formula should be replaced by S(k, e, a) 
so e appears both in S and in s. But knowledge of the complete spectral function will 
not be needed if the apparatus is capable of accepting events of a wide enough range 
of k and e so that the integral of S will be close to unity. 

Finally, we point to an important experimental factor. While it is important to 
determine k and e for each event to correct for the strong s dependence of the cross 
section, the key quantity determining the number of pure (p, 2p) events, and hence T, 
is the efficiency of the anticoincidence counters that are supposed to trigger on particles 
produced by the inelastic interactions. If these detectors remove events containing extra 
protons emitted because of the short range nucleon nucleon correlations they will lower 
T. If they allow in an event containing a pion, they will increase T, especially since 

— * 

such events will be given a high weight when k is extracted from the reconstruction. 
Thus complete knowledge of the anticoincidence efficiency is crucial to extracting T. 
This requires generating the soft collisions from the known event structures to find the 
efficiencies at the different bombarding energies. 

7 Conclusions 

We have demonstrated, using the DWI approximation, that the transparency depends 
directly on the nucleon spatial distribution p(r) and the p,p total cross section. To 
calculate T there is no need to know the transparency for each nucleon wave func- 
tion. We obtain the classical result that the p-A cross section can be factorized into 
a k independent transparency, T(er), an effective "normalized to unity" momentum 
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distribution S(k), and the p,p cross section evaluated at the appropriate Mandelstam 
variables. 

Further we point out that the DWIA calculations of Lee and Miller and our discrete 
nucleon calculation differ considerably in the magnitude of the "nuclear correlations" . 
In the latter we are dealing with finite size nucleons distributed throughout the nuclear 
volume rather than a smeared out uniform matter density, p(r). The nuclear correlation 
is incorporated simply, by not counting the struck nucleon in the absorptive path. The 
smooth correlation function used by Lee and Miller lowers the nucleon density near 
the struck nucleon in a continuous way consistent with the DWIA. But calculation 
of their C(z — z') shows that it makes for a smaller depression of p than required by 
the existence of discrete nucleons. We suggest that one should not apply the DWIA 
literally in computing T at energies of ten GeV where the proton wavelengths are 
smaller than the average internucleon spacing and the confinement of nuclear matter 
into discrete regions may better describe the physics. 

A Using a nuclear matter correlation function in 
finite nuclei. 

A nuclear matter correlation function cannot be added without further ado to a finite 
nucleus calculation, since it does not have to satisfy the condition (^). If we wish to 
remain as close in spirit to such a correlation function, we need to require a simple 
renormalization. The simplest (symmetric) way to renormalize the nuclear matter 
C NM (r) is to write 

C(n, f 2 ) = T^OC^dn - r 2 |) 7 (r 2 ). (30) 

This is symmetric under the interchange of r\ and r 2 , and we can easily solve for 7 
from the condition (0). This leads to the integral equation 

7(h) / d 3 r 2 C NM (\n - f 2 |)7(r 2 )p(f 2 )/(A - 1) = 1. (31) 

We can easily see that even if there are no correlations, C NM = 1, 7 is not equal to 1, 
(7 = y 7 ^ 1 )- This shows the finite A corrections missing in the original work of Farrar 
et al. 0. 

The integral equation (^l|) can easily be solved numerically on a grid in coordinate 
space. We have performed such calculations, and find, as expected, that the larger the 
nucleus, the smaller the effect. For 12 C 7 varies by a few percent, whereas for 208 Pb 
the effect is much smaller. 
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